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Abstract. We study the superconducting proximity effect induced in the surface states of the 
3-d topological insulator Bi2Se3 by a singlet, s-wave superconductor deposited on its surface. 
To this effect, the k • p-Hamiltonian of Bi2Se3 and the BCS-Hamiltonian are mapped onto 
tight-binding chains which we couple through a transfer-Hamiltonian at the interface. We then 
employ the Recursive Green's Function technique to obtain the local spectral function and infer 
the dispersion of the interface-states from it. In agreement with earlier microscopic studies of this 
problem, we find that the Fu-Kane model is a reasonable approximation at energies e <C Asc- 
However, for energies close to the SC bulk gap, the Fu-Kane model is expected to break down. 
Indeed, our numerical calculations show strong modifications of the interface-state dispersion 
for e > Asc- We find that the proximity effect can be strong enough to induce a gap in the 
surface state that is comparable to the superconducting gap. An analysis of the spatial profile 
of the states shows that their weight shifts towards the SC as the coupling strength increases. 
We conclude that an intermediate coupling is ideal for realising the Fu-Kane scenario. 



1. Introduction 

The research effort devoted to topological insulators has grown rapidly since their initial 
discovery in 2007 [lj and is presently developing in a multitude of directions (2J [3], one of 
them being the superconducting proximity effect. It was argued that in superconducting 
heterostructures involving 3-dimensional topological insulators, the physics of the interface 
region can effectively be described in terms of a topological superconductor, with Majorana- 
Fermions at line-junctions and vortex-cores [lj. These excitations obey non-abelian statistics 
and are thus highly relevant to the field of topological quantum computation [5J. Apart from 
proximity structures, topological superconductivity can also be realised in non-centrosymmetric 
materials, which feature spin-current carrying surface states jB]. Such surface states are 
characteristic for topological materials, where the electronic structure is characterised by non- 
trivial topological invariants. 

In their initial work, Fu and Kane described the proximity effect in a topological insulator by 
simply including a singlet, s-wave pair-potential in the effective Hamiltonian of the 2-dimensional 
surface state [lj: 





vp, A and fi can be considered as free parameters that we will use to fit the dispersion of this 
model to our numerical results. 

While the proximity effect is known to induce pairing correlations which suppress the density 
of states and, under certain circumstances, opens up a gap in the quasiparticle spectrum of 
an adjacent material, it does not induce a pair potential as such [7j. The latter would require 
a finite pairing interaction on the normal side, as can easily be seen from the gap-equation. 
Thus, there was no obvious answer to the question whether this model can be justified from a 
microscopic theory of the proximity effect, and a first investigation was provided by Stanescu et 
al. [8]. Their work was based on a tight-binding model describing the whole structure within 
a transfer-Hamiltonian approach. The authors concluded that the Fu-Kane model is a valid 
low-energy theory and that the inclusion of a pair-potential in the surface state Hamiltonian 
can be derived from integrating out the superconducting part of the structure and including it 
as a self-energy within a Green's function formalism. Since the relevant states are localised at 
the interface, this self-energy formally appears as an effective order parameter in the surface 
state Hamiltonian [9j. A similiar study has also been performed for the edge-states of 2D TIs 

A different approach was presented a little later by Lababidi et al. [11], using a numerical 
method to solve the Bogoliubov-de-Gennes equations for a TI-SC heterostructure, where the 
k • p-Hamiltonian for Bi2Se3 |X2^ IT3] was used to model the TI. For technical reasons, the 
SC-part of the structure had to be described by a two-band model, with a superconducting 
conduction band and a split-off valence band, where the conduction band is assumed to be 
identical to one of the two orbitals of the k • p-Hamiltonian. The interface is modelled by a step- 
function dependence of the Hamiltonian parameters. The authors claim that this model would 
account for the strong-coupling limit of the proximity effect, while the approach of Stanescu 
et al. was limited to weakly transmitting interfaces. They also reach the conclusion that the 
Fu-Kane model is a reasonable low-energy approximation, but do not endeavour to extract a 
surface state Hamiltonian from their calculation. They do, however, perform a self-consistent 
calculation of the SC-gap and find a suppression in proximity to the interface. 

We here propose a different approach to this problem, which can directly be applied to 
material specific k • p- or tight-binding-Hamiltonians obtained from fits to ab-initio band 
structure calculations. Thus, the starting point is similar to the work of Lababidi et al., but we 
then proceed by solving the problem with a combination of the finite differences and recursive 
Green's function methods. This has proven to be convenient tool for transport calculations in 
low-dimensional conductors |14j . By exploiting that translational invariance is only broken in 
one direction, we show that it is also a very efficient method for solving the problem discussed 
here. As opposed to the approach of Stanescu et al., this allows us to treat an infinite system, 
where potential finite size effects are absent. This is not important if the energy gaps on both 
sides of the interface are of the order of the band gap, and the interface-states are thus strongly 
localised. This assumption was made by Stanescu et al. However, a realistic SC-gap is in fact 
a lot smaller than the band-gap of the TI. The interface-states hence extend into the SC on 
the scale of the coherence length. In our model, there is also no need to maintain the same 
Hilbert-space dimensionality on both sides of the interface in the continuum models - which 
appears to be the reason for the inclusion of an additional valence band on the SC-side in the 
work of Lababidi et al. We show that the proximity induced superconducting gap in the surface 
states can in principle become comparable to the SC-bulk gap, and that the Fu-Kane model is 
a sensible low-energy theory if the coupling is not too strong. 

2. The model 

We start from the continuum model of the k • p-Hamiltonian on the Tl-side and a Bogoliubov- 
de-Gennes Hamiltonian with a single parabolic band on the SC-side. We consider an infinite 
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Figure 1. Tight-binding model of the TI-SC heterostructure. A central device, consisting of a 
one-dimensional tight-binding chain is connected to semi-infinite leads via self-energies that are 
added to the on-site terms at the respective final sites. 

system and the interface is assumed to lie in the z = 0-plane. We exploit the translational 
invariance in the x-y-directions - which implies the conservation of parallel momentum, ku, at 
the interface - by Fourier-transforming the Hamiltonians to real-space coordinates only in the 
z-direction (k z \— )■ —id z ). Using the method of finite differences [15] , we then map the model onto 
a 1-dimensional tight-binding Hamiltonian which depends parametrically on ku. The coupling 
between the two parts of the system is modelled by a transfer-Hamiltonian of the form: 

Hi = T\N + 1)(N\ + T^\N)(N + 1|, (2) 

coupling the final site ./V of the TI to the first site N + 1 of the SC via the matrix element: 

T= t - (Wti + Vsc^) . (3) 

Here, Vti,sc are the respective hopping elements of the TI and the SC derived from the 
finite differences discretisation (see Fig. [TJ. V is a projector connecting the single band of the 
SC to the two orbitals of the k • p-Hamiltonian. This allows for a greater flexibility in the 
description of orbital-coupling compared to the approach of Lababidi et al. The inclusion of a 
second 'dummy' band on the SC side is also not necessary, since we are not forced to maintain 
the same Hilbert-space dimensionality on both sides of the interface. We found that coupling 
the SC-band to the Tl-orbital with positive parity ('+'-orbital in Ref. [12J) is most effective. We 
hence consider a V here, where the matrix elements between the '+'-orbital and the SC-band 
are equal to 1 and all other elements zero. The parameter t allows us to tune the coupling 
strength. In the limit t — > 0, the two subsystems are decoupled and the surface state dispersion 
of a vacuum boundary is recovered. It agrees with the dispersion obtained directly from the 
continuum model. For t — > oo, the proximity effect is fully absent as well, but the dispersion 
is slightly modified compared to the vacuum boundary solution. To investigate the spectral 
properties of the system, we calculate the local spectral function, which is obtained with a 
recursive Green's function approach [16]. This technique allows us to describe an infinite tight- 
binding (TB) system by coupling a finite tight-binding chain to semi-infinite leads on both sides 
of the interface. The devision of the infinite TB-chain into a 'central device' and 'lead' regions 
is arbitrary, since the procedure of integrating out the leads is exact and not an approximation. 
The only requirement is that the leads are homogenous, i.e., the interface must be part of the 
central device (for the case of a homogeneous order order parameter in the SC; otherwise the 
region of varying order parameter must be included in the central region) . The Green's function 
of the central device coupled to the leads is then formally defined as: 



0(k,|, e) = (e + irj- H(k\\) - £ SC (k||, e) - E T i(k N , e))"\ 



(4) 
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Figure 2. The spectral function Ajv(kii, e) as a function of momentum and energy along the cut 
k x = k y for t = 0.92 (upper-left), t = 0.9 (upper-right), t = 0.88 (lower- left) and t = 0.86 (lower- 
right). The upper-left plot also shows a tentative fit to one branch of the Fu-Kane dispersion 
(pink line). 



where Egc and Eti are the self-energies of the SC and TI respectively, and % is the tight-binding, 
Bogoliubov-de-Gennes (BdG) Hamiltonian of the central device: 

H(k||)=% c (k||)+^ T i(k||)+^ T (k||). (5) 

The local spectral function is given by: 

M*\\> e ) = -^a lm Tr[HS(kn)|n)]. (6) 

a is the lattice spacing of the TB-model. Note that we take the full trace of Q, i.e., the Gor'kov 
Green's function, to map out the whole dispersion of the BdG-Hamiltonian. This implies the 
presence of hole- like branches in the dispersions we plot. We use a = 0.05 A and 77 = 10" 4 eV 
in the following. 

3. Results 

We use the parameters provided in Ref. [13] for the Bi2Se3 k-p-Hamiltonian, but the value of Co 
is shifted so that the chemical potential lies at the crossing point of the surface-state dispersion. 
The chemical potential of Bi2Se3 usually lies above the bottom of the conduction band, i.e., the 
material is in fact not an insulator [17} Il8] . However, chemical doping |18j or gate-control |19| 
have been reported to allow for tuning the chemical potential. For the superconductor model, 
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Figure 3. The right plot is the same as the upper-left one in Fig. [2] (t = 0.92). The left plot 
shows the same calculation with the assumptionSsc(kii, e) = Ssc(kii = 0, e = 0). The indicated 
fit to the Fu-Kane dispersion matches very well in the energy range shown here and is the same 
as the one depicted in Fig. [2j 

we assume a bulk-gap of Asc = 5 meV, a Fermi-energy E-p = 0.4 eV, and an effective mass equal 
to the free electron mass h 2 /2m = 3.81 eVA 2 . For other values of m, we obtain qualitatively the 
same results. To analyse the surface-state dispersion, we plot the spectral function A^{e, k||), 
N being the site of the TI next to the interface, in Fig. [2] for different values of the interface 
parameter t and on an energy scale comparable to the SC-gap. The band-gap of the TI is a lot 
larger than the SC-gap (about 0.3-0.35 eV [13 [18]). Spectral weight from the bulk states of the 
TI is outside the energy window shown in Fig. [2] Well-defined interface states in the TI only 
exist for energies smaller than the SC bulk-gap. Above this energy, the surface states of the 
TI hybridise with the continuum of itinerant SC-states, which results in the broadening of the 
spectral function above Asc- A tentative fit to one branch of the Fu-Kane model is indicated 
in the upper- left plot. We observe that the avoided crossing at the chemical potential of the 
interface states is shifted away from the T-point. This is consistent with the shift of fi in the 
Fu-Kane model reported by Lababidi et al. The gap induced in the dispersion of interface states 
by the superconducting proximity effect is in principle only limited by the SC-gap itself. This 
strong-coupling limit would, however, require a self-consistent calculation of the SC-gap, which 
is feasible within our approach, but we have not done it yet. Note also that the dispersion is 
strongly modified at the energy of the SC-gap and that the Fu-Kane model breaks down in this 
energy range. In the range of very strong coupling, which we reach for about t = 0.86 — 0.84 
with our choice of parameters, the interface states merge into the bulk-DOS of the SC, which 
leaks into the TI on a length scale given by the band-gap. In this case, there are basically 
no quasiparticle excitations anymore that can be considered interface states. For even smaller 
values of t the proximity effect becomes weaker again. The calculation of these dispersions is 
highly effective in our approach, as considering a single site coupled to the semi-infinite leads on 
the left and right side is sufficient to obtain A^(e, kn). When a self-consistent calculation of the 
SC order parameter is included, the TB-chain must be extended into the superconductor until 
the point where Asc recovers its bulk-value. 

While the agreement of the interface state dispersion is a necessary condition for proving 
the validity of the Fu-Kane model, it is, in principle, not a sufficient one. To complete this 
prove, an effective Hamiltonian for the interface states must be extracted from the microscopic 
calculation. This can be achieved by formulating the proximity of the superconductor in terms 
of a self-energy, i.e. coupling the final Tl-site directly to the SC-lead, and then approximating 
£gc(k||, e) as £gc(k|| = 0, e = 0). This is illustrated in Fig. |3| were we show that the dispersion 




120 



— 1=0.86 
- - -t=0.88 
— 1=0.92 



5 

z/nm 



10 



15 



Figure 4. Spatial profile of the states at the bottom of the surface-state dispersion for 
t = 0.92,0.88,0.86. Specifically, k x = k y = 0.0064 A" 1 (red), 0.013 A" 1 (green), 0.0275 A" 1 
(blue); and e = 1.13 meV (red), = 3.03meV (green), = 4.74 meV (blue). The interface is located 
at z = 3 nm. 

thus obtained agrees very well with the Fu-Kane dispersion. The off-diagonal components in 
particle-hole space of the self-energy £gc do then indeed have the structure of a singlet, s-wave 
gap. A surface-state Hamiltonian could now be obtained by considering a finite TI tight-binding 
model, were the self energy Egc(ki| = 0,e = 0) is added to the site at the interface. This TB- 
Hamiltonian can then be diagonalised exactly and projected onto the interface states, which 
would provide the interface state Hamiltonian. 

The spatial profile of the states can also be inferred from the spectral function, as shown in 
Fig. [4| For weak coupling, the interface states are mainly localised on the Tl-side, but as the 
coupling increases, their weight shifts towards the SC. The decay length on either side of the 
interface is governed by the energy Agc,Ti — e, where Ati,sc is the respective bulk energy gap, 
and e is the energy of the state. The band gap Ati is of the order O.leV, while the SC gap is 
typically two orders of magnitude smaller. Thus, we hardly see any decay on the superconducting 
side on the interatomic length-scale shown in Fig. [3j in particular when the coupling is strong 
and the band bottom of the interface state dispersion is close to the SC-gap. 

4. Conclusions 

We presented an efficient numerical approach to the proximity effect in topological surface states, 
based on the recursive Green's function technique within the transfer-Hamiltonian approach. 
The underlying tight-binding model can be derived from ab-initio bandstructure calculations 
and thus allows for material specific calculations in a straight-forward manner. We showed that 
the only fundamental limitation to the proximity induced gap in the surface states of the TI 
is the SC energy-gap itself, but an intermediate coupling regime is most desirable to achieve a 
superconducting interface state which is well localised and reasonably described by the Fu-Kane 
model at small energies. We showed that in the case of weak transmission, approximating the 
self-energy by its zero energy, zero momentum value reproduced the Fu-Kane dispersion to a 
high precision and argue that this can likely be used to extract an effective interface-Hamiltonian 
from our calculation. 
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